One-dimensional Anderson localization in certain correlated random potentials 
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We study Anderson localization of ultracold atoms in weak, one-dimensional speckle potentials, 
using perturbation theory beyond Born approximation. We show the existence of a series of sharp 
crossovers (effective mobility edges) between energy regions where localization lengths differ by 
orders of magnitude. We also point out that the correction to the Born term explicitly depends 
on the sign of the potential. Our results are in agreement with numerical calculations in a regime 
relevant for experiments. Finally, we analyze our findings in the light of a diagrammatic approach. 

PACS numbers: 03.75.-b,42.25.Dd,72.15.Rn 



I. INTRODUCTION 

Anderson localization (AL) of single electron wave 
functions [l| , first proposed to understand certain metal- 
insulator transitions, is now considered an ubiquitous 
phenomenon, which can happen for any kind of waves 
propagating in a medium with random impurities [1, 
It can be understood as a coherent interference effect of 
waves multiply scattered from random defects, yielding 
localized waves with exponential profile, and resulting 
in complete suppression of the usual diffusive transport 
associated with incoherent wave scattering So far, 
AL has been reported for light waves in diffusive me- 
dia 0,01 and photonic crystals 0,13) sound waves 0|, or 
microwaves [10|. Ultracold atoms have allowed studies 
of AL in momentum space [ill, and recently direct 
observation of localized atomic matter waves 

mm- 

In one-dimensional (ID) systems, all states are local- 
ized, and the localization length is simply proportional to 
the transport mean-free path [l3|. However, this strong 
property should not hide that long-range correlations 
can induce subtle effects in ID models of disorder, in 
particular those whose power spectrum has a finite sup- 
port [3, [l3| • Examples are random potentials resulting 
from laser speckle and used in experiments with ultracold 
atoms [ll, ill, [l^. Indeed, by construction jl^l, speckles 
have no Fourier component beyond a certain value 2fc<., 
and the Born approximation predicts no back-scattering 
and no localization for atoms with momentum hk > hk^. 
This defines an effective mobility edge at k = (l7l |. clear 
evidence of which has been reported (l3l |. 

Beyond this analysis -relevant for systems of moderate 
size [ll, study of AL in correlated potentials beyond 
the effective mobility edge requires more elaborated ap- 
proaches. In Ref. [2l|, disorder with symmetric probabil- 
ity distribution was studied, and examples were exhib- 
ited, for which exponential localization occurs even for 



k > k^ although with a much longer localization length 
than for k < k^. It was also concluded that for Gaus- 
sian disorder, there is a second effective mobility edge 
at 2k^, while for non-Gaussian disorder, it is generally 
not so. These results do not apply to speckle potentials 
whose probability distribution is asymmetric. Moreover, 
although speckle potentials are not Gaussian, they de- 
rive from the squared modulus of a Gaussian field, and, 
as we will show, the conclusions of Ref. [2l| must be re- 
examined. Hence, considering speckle potentials presents 
a twofold interest. First, they form an original class of 
non-Gaussian disorder which can inherit properties of an 
underlying Gaussian process. Second, they are easily im- 
plemented in experiments with ultracold atoms where the 
localization length can be directly measured (l3| . 




Figure 1: (Color online) Lyapunov exponent 7 calculated two 
orders beyond the Born approximation for particles in ID 
speckle potentials created with a square diffusive plate, versus 
the particle momentum hk and the strength of disorder Er = 
and ctr are the amplitude and correlation 
length of the disorder). The solid blue lines correspond to 
Er = 0.1 and Cr = 0.02. 
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In this work, we study AL in speckle potentials beyond 
the Born approximation, using perturbation theory [j^ . 
numerical calculations, and diagrammatic methods. We 
find that there exist several effective mobility edges at 
^c^'' = P^c with integer p, such that AL in the successive 
intervals fc^^"^-' < k < k^f^ results from scattering pro- 
cesses of increasing order. Effective mobility edges are 
thus characterized by sharp crossovers in the k depen- 
dence of the Lyapunov exponent (see Fig. [1]). We prove 
this for the first two effective mobility edges by calculat- 
ing the three lowest-order terms, and give general argu- 
ments for any p. In addition, we discuss the effect of odd 
terms that appear in the Born series due to the asym- 
metric probability distribution of speckle potentials. 



II. SPECKLE POTENTIALS 

Let us first recall the main properties of speckle po- 
tentials. Optical speckle is obtained by transmission of 
a laser beam through a medium with a random phase 
profile, such as a ground glass plate [131 ■ The resulting 
complex electric field f is a sum of independent random 
variables and forms a Gaussian process. In such a light 
field, atoms experience a random potential proportional 
to the intensity ■ Defining the zero of energies so that 
(y) = 0, the random potential is thus 



V{z) = V^x (|a(z/a^)p - (|a(z/a^)p)) 



(1) 



where the quantities a[u) are complex Gaussian variables 
proportional to the electric field E, and and fea- 
ture characteristic length and strength scales of the ran- 
dom potential (The precise definition of and CTr may 
depend on the model of disorder; see below). In con- 
trast, V{z) is not a Gaussian variable and its probability 
distribution is a decaying exponential, i.e. asymmetric. 
The sign of is thus relevant and can be either posi- 
tive or negative for "blue"- and "red"-detuned laser light 
respectively. However, the random potential V{z) inher- 
its properties of the underlying Gaussian field a{u). For 
instance, all potential correlators c„ are completely de- 
termined by the field correlator Ca(u) ~ {a{0)*a{u)) via 



n 



alau{i)).:{a*pau{p)), (2) 



where ap> = a{zpi /ctr) and 11 describes the p\ permuta- 



tions of {l,...,p}. Hence, C2(w) 
a{u) so that {\a{u)\'^) = 1, we have \/{V{zy) 



Ca(M)|^ and defining 



Also, since speckle results from interference between light 
waves of wavelength Al coming from a finite-size aperture 
of angular width 2a, the Fourier transform of the field 
correlator has no component beyond k^ = 27rsina/AL, 
and Ca has always a finite support: 



Ca{q) = for l^l > fc,a-R = 1. 



(3) 



As a consequence, the Fourier transform of the potential 
correlator also has a finite support: C2{q) = for |g| > 2. 



III. PHASE FORMALISM 

Consider now a particle of energy i? in a ID random 
potential V(z) with zero statistical average \V{z) need 
not be a speckle potential here] . The particle wave func- 
tion 4> can be written in phase-amplitude representation 

(j)[z) ^ r{z)s\n[6{z)]] d^cj) — kr{z) cos [9 {z)] , (4) 

which proves convenient to capture the asym ptotic decay 
of the wave function (here k = ^j2'mE/h? is the parti- 
cle wave vector in the absence of disorder). It is easily 
checked that the Schrodinger equation is then equivalent 
to the coupled equations 



d,e[z)^k[i-{v{z)/E)shi^ {e{z))\ 



(5) 



ln[r(z)/r(0)] = fc / dz' {V{z')/2E) sin (26'(z')) .(6) 

Since Eq. |[5|) is a closed equation for the phase 9, it is 
straightforward to develop the perturbation series of 9 
in increasing powers of V . Reintroducing the solutions 
at different orders into Eq. ([6]) yields the corresponding 
series for the amplitude r{z) and the Lyapunov exponent: 



7(fc) = lim 

|z|— i-CXD 



(ln[r(z)]) 



(7) 



The nth-order term thus expressed as a 

function of the n-point correlator Cn{zi, Zn-i) = 
{V{0)V{zi)...V{zn-i)) of the random potential, which we 
write C„(zi, z„_i) = Vj^c„ (zi/ctr, z„_i/crR). Up 
order n = 4, we find 



+ - / du C2(u) cos{2ku) 



where = 2ma^Vfi/h? and 

/4(«) 



(8) 



(9) 







— oo 




du / dv C3{u,v) sm{2Kv) (10) 



oo J — oo 



u rv 



du dv dw C4{u,v,w) (11) 



X {2 cos{2 Kw) + cos[2 k{v + w — u)J^ . 

Note that the compact form (fTTj) is valid provided that 
oscillating terms, which may appear from terms in C4 that 
can be factorized as C2 correlators, are appropriately reg- 
ularized at infinity. Note also that in Eq. |[8]), the coeffi- 
cients (en/fccrR)" diverge for k ^ 0, while the exact -f{k) 
remains finite for any This signals a well-known 

breakdown of the perturbative approach. Conversely, 
the perturbative expansion is valid when '-f{k) ^ k (for 
/c — > 0), i.e. when the localization length exceeds the par- 
ticle wavelength, a physically satisfactory criterion. 
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Figure 2: (Color online) Functions /„ for ti = 2, 3 and 4 for 
a speckle potential created with a square diffusive plate [solid 
lines; see Eqs. lfT2|l and l[T3|l and the Appendix] and compari- 
son with numerical calculations (points with error bars). The 
inset is a magnification of function /4 around k = 2. 
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Figure 3: (Color online) Lyapunov exponent 7(fc) versus the 
particle momentum k as determined numerically (solid red 
lines) and by perturbation theory up to order 4 (dashed blue 
lines), for a speckle potential created with a square plate. The 
dotted green lines are the Born term. Inset: comparison of 
odd and even contributions in the Born series for Cr = 0.1. 



IV. ONE-DIMENSIONAL ANDERSON 
LOCALIZATION IN SPECKLE POTENTIALS 

A. Analytic results 

Let us now examine the consequences of the peculiar 
properties of speckle potentials in the light of the above 
perturbative approach. For clarity, we restrict ourselves 
to ID speckle potentials created by square diffusive plates 
as in Refs. for which Ca{u) = sin(u)/ti and Ca(q) oc 

0(1 — \q\) where 6 is the Heaviside step function [24 |. 
Using Eqs. ^ and ((lOl), we find 

f2{K) = ^e(i-^)(i-«;) (12) 

Mk) = -^e{l - K)[{l-K)\n{l-K) + K\n(K)] (13) 

The functions /2 and /s are simple and vanish for k > 1 
(see Fig. [2]). This property is responsible for the existence 
of the first effective mobility edge at k = such 
that j{k)a^ i^n/ka-fi)'^ for k < cr~^ while 7(fc)(TR = 
0(eR/fc(Tji)* for k > (t~^ . The fact that /a vanishes in 
the same interval (k > 1) as /2 exemplifies the general 
property that odd-n terms cannot be leading terms in 
any range of k because j{k) must be positive whatever 
the sign of V^. For k < 1 however, fsin) is not identically 
zero owing to the asymmetric probability distribution in 
speckle potentials. The term 7^^^ can thus be either pos- 
itive or negative depending on the sign of [22| • 

The function is found similarly from Eq. ifTTj) . While 
its expression is quite complicated (see the Appendix), 
its behavior is clear when plotted (see Fig. [21). Let us 
emphasize some of its important features. First, there 
is a discontinuity of the derivative of /4 at k = 1/2. 
Second, we find a very narrow logarithmic divergence, 
/4(k) ^ — (7r/32)ln|l — k| at k 1, which signals a 
singularity of the perturbative approach (note that it 



does not appear in Fig. [T] due to finite resolution of 
the plot). Finally, the value k = 2 corresponds to the 
boundary of the support of /4, showing explicitly the ex- 
istence of a second effective mobility edge at fc = 2a~^. 
Hence, while 7(A:)(Tp{ ~ (ea/ZccTii)^ for k < 2a~^, we have 
"f{k)<7^ = 0{e^/k(j^)^ for k > 2a~^, since /4(k) as well 
as /5(k) vanish for k>2. 



B. Numerics 

In order to test the vaHdity of the perturbative ap- 
proach for experimentally relevant parameters, we have 
performed numerical calculations using a transfer matrix 
approach. The results are plotted in Fig. [H Cr = 0.02 
corresponds to V^/h = 2tt x 16Hz in Fig. 3 of Ref. 
and en = 0.1 to V^/h = 27r x 80Hz in Fig. 3 and to 
Fig. 4 of Ref. [l^. For = 0.02, the agreement between 
analytical and numerical results is excellent. The effec- 
tive mobility edge at fc = cr"^ is very clear: we find a 
sharp step for j{k) of about 2 orders of magnitude. For 

= 0.1, we find the same trend but with a smoother and 
smaller step (about one order of magnitude) . In this case, 
although the Born term for k < and the fourth-order 
term for k > provide reasonable estimates (within a 
factor of 2), higher-order terms -which may depend on 
the sign of V^^- contribute significantly. 

The contribution of the odd terms can be extracted by 
taking 7+ — 7", where 7^ are the Lyapunov exponents 
obtained for positive and negative disorder amplitude of 
same modulus jVal, respectively. As shown in the inset 
of Fig. [31 the odd terms range from 30% to 70% of the 
Born term for 0.6 < ka^ < 0.9 and Er = 0.1, and are of 
the order of 7^^^ in weak disorder and away from the di- 
vergence at k ~ . This shows that the first correction 
7*-'^-' to the Born term can be relevant in experiments. 

For completeness, we have calculated the /«(«;) as the 
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S^- ^^-^ (». <» f; 

a) * b) ^ — ^ c) d) ^ ^ 

Figure 4: Relevant fourth-order backscattering contributions. 
Contrary to the case of uncorrelated potentials (26l. [27l|. the 
sum of diagrams (a)-(c) does not give zero for speckle poten- 
tials; only diagrams (b) and (d) contribute for ka^ G [1, 2]. 

coefficients of fits in powers of tT^/ka^ using series of cal- 
culations of 7(fc) at fixed k and various Er. As shown in 
Fig. [21 the agreement with the analytic formulas is excel- 
lent. In particular, the numerics reproduce the predicted 
kink at k = 1/2. The logarithmic singularity around 
K = 1 being very narrow, we did not attempt to study it. 



V. DIAGRAMMATIC ANALYSIS 

Let us finally complete our analysis using diagram- 
matic methods, which allow us to exhibit momentum ex- 
change in scattering processes as compact graphics, and 
thus to identify effective mobility edges in a quite gen- 
eral way. In ID, the localization length can be calculated 
from the backscattering probability of using quan- 

tum transport theory. The irreducible diagrams of ele- 
mentary scattering processes in speckle potentials have 
been identified in Ref. jl^. 

To lowest order in (Born approximation), the aver- 
age intensity of a plane wave with wave vector k backscat- 
tered by the random potential is described by 

fe ->»^ k k ->«K- k 

U2{k)^q + k~i; 'f.q - k =: r2k. (14) 

fc k k k 

The upper part of the diagram represents ip (particle) and 

the lower part its conjugate ip* (hole). The dotted line 

q 

•■■■■>■■■■• = eRCa(g) represents the field correlator; simple 
closed loops over field correlations can be written as a 

potential correlation » 8. Backscattering requires 

diagram ([14)) to channel a momentum 2k, entering at 
the particle, down along the potential correlations to the 
hole. Therefore, the diagram vanishes for ka^^ > 1. 
At order e^, the only possible contribution is 

.•■<■•. 

k -»e >— *<- k 

■• P :' 

Us^k) = q + k\ ^q-k + C.C. (15) 

The straight black line stands for the particle propaga- 
tor [Ek — Ep + iO]~^ at intermediate momentum p. Dia- 
gram l|15p features two vertical field correlation lines, just 
as diagram ([14)) , and thus vanishes at the same threshold 
k ~ <T~^. Evaluating two- loop diagram ([15)) . we recover 
precisely contribution ([T3)) . 



Many diagrams contribute to order e^. First there are 
the usual backscattering contributions with pure inten- 
sity correlations [Fig.[4i;a)-[4l[c)]. Both Figs.[4i;a) and[4i;c) 
have a single vertical intensity correlation and vanish for 
k > a~^. In contrast, the crossed diagram [DJb)] has 
two vertical intensity correlation lines and can thus ac- 
commodate momenta up to fc = 2a~^. Performing the 
integration, we find that this diagram reproduces those 
contributions to /^{k) for n £ [1,2] that contain factor- 
ized correlators (see the Appendix). Second, there are 
nine more diagrams, all with non-factorizable field corre- 
lations [iHl- A single one has not two, but four vertical 
field correlation lines, shown in Fig.|4ljd), and contributes 
for fco-R e [lj2]. Carrying out the three-loop integration, 
we recover exactly the non-factorizable contributions to 
for Kg [1,2]. 



VI. CONCLUSION 

We have developed perturbative and diagrammatic ap- 
proaches beyond the Born approximation, suitable to 
study ID AL in correlated disorder with possibly asym- 
metric probability distribution. In speckle potentials, the 
k dependence of the Lyapunov exponent exhibits sharp 
crossovers (effective mobility edges) separating regions 
where AL is due to scattering processes of increasing 
order. We have shown it explicitly for k — and 
k = 2ct~^, and we infer that there is a series of effective 
mobility edges at fc = pa~^ with integer p since, gener- 
ically, diagrams with 2p field correlations or p intensity 
correlations can contribute up to fc = pcr~^ . This is be- 
cause, although speckles are not Gaussian, they derive 
from a Gaussian field. Finally, exact numerics support 
our analysis for experimentally relevant parameters, and 
indicate the necessity to use higher-order terms in the 
Born series, even for fc < cr^^. Hence, important fea- 
tures that we have pointed out, such as odd terms in the 
Born series for fc < cr~^ and exponential localization for 
fc > cr~^, should be observable experimentally. 
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APPENDIX 

Here, we give the explicit formula of the function fi{K) 
for a speckle potential created by a square diffusive 
plate, such that the fourth-order term in the Born 
expansion of the Lyapunov exponent 7 reads 7^^-' = 

( to") fii^'^R)- The function fiin) is the sum of 
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three terms with different supports, where f " {k) lives on the interval k G [a,f3], and 



/f'il(K) = ^ |4-6k- -2k) - (4-2K)ln(K) - (^^ - ln(l - k) + (^^ + + k) 

-(4 - 8k) In^(K) + 22(1 - k) ln^(l - n) + (18 + 14k) ln^(l + k) 
-16(1 - k) ln(l - k) ln(K) - 4(1 - k) ln(l - k) ln(l + k) - 32(1 + k) ln(K) ln(l + k) 

-24(1 + K)Li2(K) + 32(1 + K)Li2 ( — ^ I -SkLIs ( | - 8(1 - 2K)Li2 [2-- 

\1 + kJ \1 + kJ \ K 

/f '^'(k) = ^|"^+(^ + y) '^' + 4KLi2(l - k) - f--2 + K\ ln(K - 1) - 2(k- l)ln2(K - 1) + 4Kln(K- 1) ln(K) 



r 



where Li2(z) = j" dt ln(l - t)/t = E^i -^V^^ is the 
dilogarithm function. 
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